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Abstract: We study the lattice effects on the simple holographic toy model; massive 
U{1) gauge theory for the bulk action. The mass term is for the U{1) gauge symmetry 
breaking in the bulk. Without the lattice, the AC conductivity of this model shows 
similar results to the holographic superconductor with the energy gap. On this model, 
we introduce the lattice effects, which induce the periodic potential and break the 
translational invariance of the boundary field theory. Without the lattice, due to 
the translational invariance and the mass term, there is a delta function peak at 
zero frequency on the AC conductivity. We study how this delta function peak is 
influenced by the lattice effects, which we introduce perturbatively. In the probe 
limit, we evaluate the perturbative corrections to the conductivities at very small 
frequency limit. We find that the delta function peak remains, even after the lattice 
effects are introduced, although its weight reduces perturbatively. We also study the 
lattice wavenumber dependence of this weight. Our result suggests that in the U{1) 
symmetry breaking phase, the delta function peak is stable against the lattice effects 
at least perturbatively. 
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1. Introduction 

Finding the high Tc superconductors on the cuprates system, which are not de- 
scribed by the usual BCS theory, is a remarkable breakthrough in the condensed 
matter physics occurred almost 26 years ago. Even though there are huge devel- 
opments after that, both theoretically and experimentally, we are still missing the 
core mechanism governing the system. One of the common mysterious phenomena 
in these superconductivity is its Non-Fermi-liquid behavior above the superconduct- 
ing phase, which occurs near the quantum critical point. The main difficulty to 
understand the mechanism of these high superconductor is, of course, its strongly 
coupled dynamics and its lack of the normal quasi-particle pictures. We do not 
yet fully understand by what mechanism and for what materials, how high Tc su- 
perconductor can occur in nature. However, as is recent discovery of iron-based 
superconductor, experimental progresses on this field are remarkable. These include 
recent developments of the cold atom experiments. Therefore it can happen that in 
the near future, we get more crucial experimental data which helps to deepen our 
understanding of the core mechanism. 

On the other hand, one of the most surprising development coming from the 
string theory is the realization of holographic principle, which states that two totally 
different theories, string (or gravitational) theories in asymptotically anti-de Sitter 
space background and strongly coupled large gauge theories, are equivalent at 
some limit [1, 2, 3]. Recent developments of the holography by applying that to 
the strongly coupled condensed matter system, is just tremendous^. Especially the 

^See for examples, [4, 5, 6, 7] 
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construction of the holographic superconductor (superfluid) [8, 9, 10], where the 
U{1) symmetry breaking through the hairy black hole in the bulk, intrigues the 
many interesting developments. See, for examples, [11, 12] for the review of the 
holographic superconductors. 

In the real-world materials, it happens quite frequently that the materials show- 
ing the superconducting phase do not have a translational invariance and Lorentz 
symmetry, due to the crystal structure of the background atoms. In high super- 
conductor, the effects of the background atomic structure are very important and it 
is expected that two-dimensional structure plays the significant role. One of the very 
important effects of the background atomic lattice is that it violates the translational 
invariance and Lorentz invariance and induces the periodic potential. 

If the material possesses a translational invariance and a net charge, it is more 
or less guaranteed that its electric conductivity shows the delta function peak at 
zero frequency. This is simply because of the fact that the charged objects are 
kept accelerated by the outer electric field in a translationally invariant system. 
We can also understand this from the fact that by the Lorentz boost, the system 
acquires a nonzero current with zero applied electric field. However, once we break 
the translational invariance by the lattice, then, there is no guarantee that such a 
delta function peak appears on the conductivity^. Therefore it is quite interesting 
to consider how the delta function peak in many interesting holographic system is 
influenced by these lattice effects. In this paper, we take a first step towards the 
lattice effects on the holographic superconductor; we study the lattice effects on a 
toy model, which is massive f/(l) gauge theory for the bulk action. The mass term 
of gauge boson is for the U{1) gauge symmetry breaking in the bulk. 

Our toy model, although it is a different theory from the holographic supercon- 
ductor model analyzed by [9, 10], has properties which is similar to the holographic 
superconductor; Without the lattice, it shows the mass gap, and AC conductivity 
is quite similar to the results of [9]. It also shows the delta function peak at zero 
frequency. Therefore, we find it interesting to ask how the zero frequency delta func- 
tion peak in our model, which is due to the translational invariance and the mass 
term, is influenced by the lattice effects. Since in our model the gauge boson has 
mass term, corresponding to the U{1) symmetry breaking of the superconductor (su- 
perfluid) phase, this analysis is a first step to study the generic lattice effects on the 
holographic superconductor phase. 

There are several technical points which are worth quoted at this stage. In this 
paper, we consider the probe limit, namely we neglect the effects of the gravity. It 
is known that in the normal phase {i.e., non-superconducting phase) without taking 

^Even if the system has a translational invariance, if we apply magnetic field, this delta function 
peak disappears. This is because momentum is not conserved in the presence of magnetic field. 
The same is true in the presence of disorder. See, for example, [13, 14] for the explicit examples of 
these in the holographic conductivity calculations. 
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into account the gravity, conductivity becomes trivial and there is no delta function 
peak appearing. The delta function peak appears only if we take into account the 
gravity effects in the normal phase. On the other hand, in the superconducting phase 
of [9] , the delta function peak appears without taking into account the gravity effects, 
so is our bulk massive U (1) gauge model. In this paper, without taking into account 
the gravity, we study how the delta function peak in our model is influenced if we 
introduce the lattice effects perturbatively (periodic disturbance to the system). 

Before we end this introduction, we comment on several closely related refer- 
ences. Recently there are developments to take into account the lattice effects for 
the holographic condensed matter system. In the paper by Maeda, Okamura, and 
Koga [15], the geometry, where the back reaction of lattice effects is taken into ac- 
count perturbatively, is constructed. In that paper, the lattice effects are introduced 
through the chemical potential. In [16], Horowitz, Santos, and Tong calculated the 
conductivity under the presence of the lattice for the normal phase, namely, non- 
superconducting phase. They showed, using a very powerful numerical technique, 
how the AC conductivity, especially the zero frequency delta function peak, is influ- 
enced by the lattice effects. They showed that the delta function peak disappears 
by the lattice effects, as is expected for the properties of the normal phase. In that 
paper, the lattice effects are introduced through the neutral scalar field. Very re- 
cently, a paper [17] by Liu, Schalm, Sun, Zaanen, appears where on the geometry 
without gravitational back reaction, they discussed the lattice effects for the holo- 
graphic fermion correlators, especially its pole for the Non-fermi-liquids, to see how 
their dissipation relations are modified. 

2. Massive U{1) gauge boson model 
2.1 The model 

In this paper we consider the following toy model of holographic superconductor for 
the bulk action. 



where V{u, x) is external potential, and plays the role of position-dependent mass 
term for the gauge boson A^. u is radial coordinate in the bulk and x is the spatial 
coordinate in the boundary theory. In more realistic holographic superconductor 
model, V{u, x) is given by the condensation of the charged scalar field as [8, 9, 10]. 
There, the massless U{1) gauge boson couples to the charged scalar \E', where takes 

non-zero VEV, -qf^ackground 
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This gives the potential 

l^(u,a;) ~ 1^^""'^'?™''"^!^ (2.3) 

corresponding to the spontaneous U{1) symmetry breaking in the bulk, which is 
dual to the global U{1) symmetry breaking in the boundary theory. However in this 
paper, we consider V{u, x) as given input for the symmetry breaking. Especially, we 
consider V{u,x) which satisfies 

V{u,x) -> (at the boundary), (2.4) 

so that at the boundary, the mass term V{u, x) for the gauge boson disappears. We 
would like to calculate the AC conductivity through the bulk f/(l) dynamics A^, and 
discuss the delta function peak with the lattice effects. 

One of the main reasons why we consider this toy model is its simplicity. By 
restricting the degrees of freedom, the calculation for the AC conductivity, especially 
by choosing appropriate boundary condition, becomes much simpler than the holo- 
graphic superconductor model in [9, 10]. However, as we will show later, this model, 
without the lattice, shows the AC conductivity which is quite similar to the one of 
the holographic superconductor model. It shows the energy gap. Since one of the 
essential features of the holographic superconductor model is the U{1) symmetry 
breaking in the bulk, we expect this bottom-up model captures some of the essential 
features. In this paper, we consider the lattice effects on this model. 

The equations of motion for gauge field become 

V.F'^^ = -^d,{,/^Fn = -2^(u, x)A^' , F^, = V,A, - V,A, . (2.5) 

We give, by hand, the non-zero arbitrary VEV for the background gauge potential 
Af^ as 

A^ = x) ^ , (2.6) 

which should be correlated with the nonzero VEV of the charged scalar field (2.2) in 
holographic superconductor [9, 10]. This nonzero VEV produces a net charge for the 
boundary theory through the normalizable mode of j\^^(^^9-round ^ q-^ background, 
in order to calculate the conductivity, we will add small fiuctuations given by the 
ansatz 

SA^dx^" = SAt{t, u, x)dt + M^.(t, u, x)dx + M„(t, u, x)du . (2.7) 

From the dynamics of these fiuctuations 5A^, we would like to calculate the AC 
conductivity. However, since the U{1) gauge boson equations of motion are linear 
equations, all of our analysis are independent on the VEV of the gauge boson (2.6). 



-4 - 



Therefore our argument for the calculations of the conductivity is independent on 
the background VEV (2.6). 

We take the background metric to be 



corresponding to the Schwarzschild AdS black brane. 

We study the lattice effects on our model given by the action (2.1). As we quoted 
already, this model is a different theory from the holographic superconductor one [9] , 
where the degrees of freedom is U{1) gauge boson and the charged scalar field. We 
give the U{1) symmetry breaking background by hand as (2.2) and (2.6). We have 
not taken into account the dynamics of the charged scalar field fluctuation in the 
f/(l) symmetry breaking phase, which might not be consistent from the original 
gauge boson and charged scalar field theory in the bulk. Therefore the reader may 
regard our model a rather toy model akin to the U{1) symmetry breaking phase of 
the holographic superconductor. Although our massive U{1) gauge boson model is 
different model, there are several merit to study this model. First, due to the lack 
of the dynamical charged scalar, the analysis of solving the equations of motion is 
simpler in this model. Second, even though we will not take into account the charged 
scalar field fluctuation for the conductivity calculation, the results of conductivity 
without the lattice in our model are quite similar to the superconductor model, it 
shows the mass gap. Technically, this is because charged scalar fluctuation does not 
directly couple to the gauge boson fluctuation without the lattice. Third, this simple 
model also shows the zero frequency delta function peak in the probe limit of the 
holographic superconductor. Therefore even though this model is different, since 
this model shows very similar properties to the U{1) symmetry breaking phase of 
the holographic superconductor, we regard this model belongs to the same category 
to the one in [9]. Therefore even in this model, there are interesting questions we can 
ask. In this paper, we restrict our attention to this model, and we study the lattice 
effects on it. 

We also point out that in this work, we neglect gravitational back reaction. We 
expect that even if we go beyond the probe limit, the back reaction does not change 
the picture drastically as the case of holographic superconductivity. For this reason, 
we take the background to be Schwarzschild AdS black brane as (2.8), even though 
there are non-trivial background flux (2.6). Once we take into the back reaction, the 
metric will be modified either into Reissner Nordstrom AdS type of brane or AdS 
hairy black brane solution. 




(2.8) 



where h is radial dependent function, given by 



h{u) = 1 - , 



(2.9) 
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2.2 Generic analysis 

We would like to consider the gauge boson fluctuation so that we can obtain AC 
conductivity. For each of the fluctuation (2.5), the equations of motion becomes 

hSAt,uu - {SA.^tx - SAt^,,) - h 6A^^^t -V{u,x)6At = , (2.10) 

8 8 2/^/? 

h—{h6A^,u) - {6A^,u-SAt,,t)-h—{h6Au,x) -Viu,x)6A, = 0, (2.11) 

8u 8u 

2L'^h 

5At,tu - hdA^^^u - 5Au,tt + h5Au,xx :^V{u,x)5Au = . (2.12) 

Equation of motion for the y component is trivially satisfled by 5 Ay = 0. 

Clearly, if the potential V{u, x) is independent on position x, then there exists 
the solution where both 6 At = 6Au = with nontrivial SA^, which is independent 
on X. In such cases, fluctuation equation for SA^ becomes single differential equation 
for the second order. 

Taking the simple time-dependence as 6Ai = e~^'^^ai {i = t, x, u) for the AC 
conductivity, we obtain 

2L2 

hat^uu + at.xx + + ^w/ia„,„ ^y{u, x)at = , (2.13) 

8 8 2L?'h 

h^{hax,u) + u}^ax - iujat,x - h—{hau,x) ^^iu, x)a^ = , (2.14) 

8u 8u 

2L'^h 

- iuat,u - hax,xu + co'^au + hau,xx ^Vi^, x)au = . (2.15) 

We would like to solve these coupled differential equations by perturbation and 
obtain the AC conductivity. For that purpose, we add the lattice effects which simply 
take the following cosine form by the perturbation as 

V=^ {Vo{u) + e6V{u) cosgx) , (2.16) 

where e is a small parameter. Vq{u) corresponds to the homogeneous charged scalar 
condensation, while 6V{u) corresponds to the lattice effects. We take the ansatz that 
the lattice has position x dependence given by the wavenumber q. We will take an 
explicit example later, but for a moment we keep it generic w-dependent functions 
Vo(m) and 6V{u). 

We will conduct perturbation expansion for small e as 



ttx = a^x^ (w, u) + ea^^^ (w, u) cos qx + e^a^^-* (w, m, x) + • • • , (2-17) 
at = ea[^^ (w, u) sin qx + e^af* (w, u, x) + ■ ■ ■ , (2. 18) 

ttu = ea^^-* (w, u) sin qx + e^a^^^ (w, m, x) + ■ ■ ■ . (2. 19) 
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Note that af'^ and a[^^ are independent on x but we keep the imphcit x dependence 
for af\ Then, from (2.13), (2.14), (2.15), we obtain 

ft|^a,"> - (,'ai'> + i,a,a«) + ii^ft|^a,<'l - '^a? = 0, (2,21) 

- >^|;a!" + ,A Aai" + ^=a<') - g'h a« - H^aW = , (2.23) 

, d"^ (2) , (2) , . "9 (2) . ■ , d 2\4) (2) 25\^ (i) . 

+ + ^^^"'x + ^^h—a\ ^ -a\ = -^al sm qx cos qx , 

oti^ ox^ ox ou 

(2.24) 

, d ( , daf^\ 2 (2) -9 (2) , d / , 9al^M 2hVo (n) 2 h6V cos"^ qx 
'a^I 1^'^ j + - 4 ^ - -^a^ - [h^ j - -^ai = 4 , 

(2.25) 

d (2) &^ (2) 2 {2\ 2/iVo (2) 2h6V /-i) 

- luj—al - h——a\. ' + h—^a\, ' + uj a\,' ^a}, ' = — ^a}, ' sm qx cos qx . 

du dudx dx^ 

(2.26) 



In order to simphfy, we define the average physical quantities over the spatial 
direction x on the range 2n/q as 

A{u) := — / A{u,x)dx. (2.27) 
271" Jo 

Then, due to the periodicity of the perturbation along the x direction, (2.25) becomes 
a simple differential equation as 

d ( ^da^x^\ 2~{2) 'ihVo^) h5V 
^Tu j + - «^ ^ - =^^^'^ ^2.28) 

namely, ai'' decouples from af^, au\ 

Furthermore in order to impose ingoing boundary condition at the horizon, we 
re-define the fields as 

^(n) _ ^^uu,^in)^ n = 0, 1, 2, ■ ■ • (2.29) 

for i = (t,x,u), where 

r du , , 
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Then, from (2.20), (2.21), (2.22), (2.28), we obtain 

.1^^) + + 2..) ^f-^ef^O, (2.31) 



dw^ au 



icuh' Lu 2 2Vc 

q 



h h ^ 

,(1) 



St 



ig^ei') + ..^{W - i^A^ , (2.32) 



du^ du 

.gc + g/.') + qh'-^ + ^d^^ + (2.33) 



yd" - + + + - i-h - «i" = , (2.34) 

h'^ + (h- + 2,u)^-?^W^=%e. (2.35) 

where ' means the w-derivatives. The zeroth order ^j"'' gives the conductivity without 
the lattice effects. By solving (2.31), we can obtain the zeroth order conductivity, 
i.e., conductivity without the lattice. 

Let us first concentrate on the leading perturbation We are interested in 
the zero frequency delta function peak. For that purpose, it is enough to study the 
behavior of these equations at low frequency limit. For that purpose, we expand 
those fields in small u as 

ef^ = ef)'° + c.eP + a;^eP + --- (2.36) 

for i = t,x, u. 

Then at the 0{uj^) order corresponding to the static limit, we have 

hf-,&' + ft'^ - ^^i"'" = 9ft'fi"^° + 9/!^ + (2.38) 



qh 



du^^"" du U2 ^ ^" ' ^ du U2 ' 



|^+(-,^.-^)d-^0. (2.39) 

Therefore in this static limit a; — > 0, ^t^^'" decouples from ^i^'''" and ^u^''^ . Then it 
is determined by solving the equation (2.37), which gives unique radial coordinate 
u dependent solution once we give the following two boundary condition; The non- 
normalizable mode of must vanish at the boundary. must vanish at the 

horizon, so that Wilson loop 

e/^M^x" _ g/A.dr ^2.40) 
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vanishes on the trivial cycle at the horizon, in the Euclid signature. These boundary 
condition determines that 



= . (2.41) 



Given this, for the next order 0{u^), we have 



dw^ du du 

= + gh'^^'^ + qh^ + '\ (2.43) 

From (2.39) and (2.44), we can write down ^u^''^, ^u^'^ in terms of ^^^'^ as 

jf(l).0 

= 1 (2 45) 

f = ( + i£(i).o ) (2 46) 

^" -2 , 2M, + }^^^ I • 



„2 



These give 



From this, it is straightforward to check that at the horizon where /i — )■ 0, ^^^^ is 
divergent-free.^ 

By plugging this (2.47) back to the equations (2.32) and (2.33), we obtain dif- 
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To see this, note that by using the regularity of £,t at the horizon, (2.42) gives 

(2.48) 



So we have 



e' = --^^6'''" + 0{.'). (2.49) 
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ferential equations for and up to order 0{u ) accuracy as 



2 ^^0 

-q ^ 



dv? 



(2.53) 



^ q {iquj + qh') ( ] , l A ( ? ( , 

+ ^ rf^ + Z,^- ) ^^^du \q^+'-^ \ du ^ h^^ 

+ + 0(c.^), (2.54) 

Note that terms like OJ^^^'^ are 0(c<;^). We can solve those equations and resultantly 
we can determine the conductivity in the low frequency limit. The results allow us 
to check if zero frequency delta function peak exists or not at this order. 
Conductivity is given by 

. (2.55) 

which is generically position dependent. Here in numerator, dominant term is a 
normalizable mode, and in the denominator, dominant term is non-normalizable 
mode. Expanding and At by e, we obtain. 



a = a^^\uj) + e cos qxa^^\uj) + e^a^'^\uj,x) + ■■■ (2.56) 



Therefore, we have, at the horizon where ft, — > 0, neglecting O(w'), 

Here, by using the similar argument, from (2.34), we can see diverges at the horizon as 

eW^'=0(ft-i). (2.51) 

Let's consider the leading divergent terms in above £,ij.S,u9^^ ■ Under the regularity condition for ^2;, 
the leading divergent terms which blow up as 0{h~^) arc 

u.^^- - "^.h-w ie^^r + (^^^e^-") + oih^) . (2.52) 

So at least up to O(aj^), the leading divergent terms, which behave as in ^^^1^5^", cancel. 
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and 



= 1 + . (2-57) 

where ' is w-derivative. If we look at the spatially averaged part of the conductivity 
a as 



adx 

= + ea(a;)(^) + e^a{ujf^'^ + • • • , (2.59) 

then we obtain, 

a^^\u) = 0, 



-(2). ^ ^ (0) _ &mT{^) , q^i'\o)6'\o) , ^^^\or6'^' jo) 

^ ^ zu;^i'\0) 2tu;e'm' 2tu;e'KQr 

^q_i!iMM _ g^ef^(0)-ei°^'(0) .2 60) 

2a;2ef (0)2 2ia;3ef (0)^^ 



We would like to evaluate these perturbative corrections to the conductivity by the 
lattice effects. However, in the real-world experiments, we usually apply a homo- 
geneous electric field and see the conductivity. This means, that we should choose 
the boundary condition such that inhomogeneous parts of the electric field are set 
to be zero. This corresponds to choosing non-normalizable modes of the 0(e) terms, 
^i^\u) for (i = t,x), are set to zero 

ef^(0)=0 , {i^t,x), (2.61) 

since they have cos qx dependence. Therefore, above conductivity formula reduces 
to 



These are the quantities which we will evaluate. 

Without solving the equations of motion explicitly, we can guess how the solution 
behaves at the zero frequency limit, and therefore, how the zero frequency delta 
function peak behaves at a; — > limit at this stage. 
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In the — )■ limit, as we showed in (2.36) and (2.41), we have 

d'^ = 0(u;) , ef = 0(l) , e^ = 0{l). (2.63) 
Similarly we can also confirm that 

= 0(1) . (2.64) 
Then by plugging these into the perturbative results (2.62), we can obtain that 

Inia^o) ~ - , Ima(^) ~ - , Ima^^^ ~ - . (2.65) 

U! U CO 

Therefore we can guess that imaginary of a has a simple pole structure ~ ^, as 
far as we consider the lattice effect perturbatively. However this argument is very 
naive, since there is a possibility that the residue of the pole becomes zero and pole 
disappears. For example, in the normal phase without lattice structure, the above 
argument breaks down since becomes zero therefore the residue of the ^ vanish. 
Therefore in order to confirm above expectation, we will now solve the equations of 
motion in more explicitly and obtain the solutions. Then, we would like to see if the 
imaginary of a has a pole as 

Im(T~-. (2.66) 
Once this behavior is confirmed, using the Kramers-Kronig relation 

H.Ml = -ipr<i.'5!M^, (2.67) 

TT J^^ Uj' -UJ 

we can conclude that there is a zero frequency delta function peak. However even if 
the delta function peak remains, how its residue (weight), u x Ima changes as we 
vary q is very nontrivial. We will see through the numerical analysis that the weight, 
which is X Ima, decreases as we increase q. 

2.3 Explicit examples for AC conductivities at low frequency limit 

In order to work on some explicit examples, let us consider examples where we take 

Vo = diu^l + d2U^) , 6V = u\ (2.68) 

Now we numerically calculate the conductivity a^^\uj), a^^\uj), a^'^\uj), given by the 
formula (2.57) and (2.62) for various q in the two typical cases, (i) di = 5, d2 = 
and (ii) di = 5, d2 = 2, respectively. We have chosen the power for Vq as and u^. 
These choices are due to the fact that the asymptotic behaviors of the background 
charged scalar, which condensates, behaves as ■qjba^'^ground ^ ^ or ^t'^o-^^a-round ^ 

see eq. (8) of [9]. 
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Figure 1: Im[cr^'^^] is plotted for various u (red) in the case (i). The best fitting curve (blue) 
is Im[(T(o)(w)] ~ -9.4 x 10"^ + 3.1/u. 



By imposing regularity condition at the horizon, 

^i°Hl) = regular, (2.69) 

we obtain the conductivity a^^\uj) by numerically solving (2.31). Fig. 1-4 show 
a^^\uj). As is seen from the figures for the Re[(T'^°''(ci;)], for the both cases, the gap 
appears. The energy gap in Fig. 4 is larger than the one in Fig. 2 because the 
"condensation" Vq{u) ~ ^qjbackground^2 -j-^ ^YiQ case (ii) is larger than the one in the 
case (i). As we mentioned in the introduction, these results for conductivities are 
quite similar to the conductivity calculations for the holographic superconductor [9] . 
Especially there are energy gaps and delta function peak at Re[o"'''^^(a;)] = 0, which is 
seen from the imaginary parts of ctq behave as K For in Fig. 2 and 4, both Re[cr^°^(ci;)] 
approaches small but nonzero value at w — 0, corresponding to the existence of the 
mass gap. 

By Kramers-Kronig relation (2.67), the real part of the conductivity contains 
a delta function peak such as Re[(T(°^(a;)] ~ 7rC*^°^5(c<;) if the imaginary part of the 
conductivity contains a pole such as Im[(T(°)(ci;)] ~ C^^yu. Fig. 1 and 3 show that 
the pole exists at a; = for both cases (i) and (ii). The best fitting curves determine 
the coefficients (7^°) as C(o) = 3.1 (the case (i)) and = 3.25 (the case (ii)). 

The first order conductivity a'^^\uj) is obtained by numerically solving the fol- 
lowing equations from (2.53) and (2.54); 



(1) iuhq (P'ix^^ 4iujhqdid2U ci^i^^ ^ 2^ 



'^"^^^ ~ g2 + 2di(l + 4n2) du^ {q^ + 2di{l + d2U^)y du 



(2.70) 
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Recr, 




Figure 2: Re[cT('''*] is plotted for various a; in the case (i). There is a delta function peak 
at a; = 0. 
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Figure 3: Im[(T('^)] is plotted for various uj (red) in the case (ii) 
curve (blue) is Im[o-(o)(w)] ~ -8.4 x 10"^ + 3.25/a;. 



The best fitting 



RecTr 



Figure 4: Re[(T^'^-'] is plotted for various uj in the case (ii). There is a delta function peak 
at (J = 0. 
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(2.71) 
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Figure 5: (color online) Im[(T((x>)(^)] in the case (i) for q = 2. The best fitting curve is 
o-(i) ~9.2 X 10^^ + 0.41/a; 
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Figure 6: (color online) Im[(7(a;)(^''] in the case (i) for q = 5. The best fitting curve is 
o-(i) ~ 0.0024 + 0.78/lj 



To obtain the solution, we need to impose boundary conditions both at the 
horizon and the infinity. As we mentioned, in real-world experiments, we usually 
apply a homogeneous electric field and measure the conductivity. Therefore we shall 
impose a constant electric field condition at infinity. Since 0(e) part of the flux has 
cosgx dependence, we require that 



Ei^\0) = (iuja^^\0) + q4^\0)) cosqx = , 



(2.72) 



for arbitrary x. We have chosen ^i^^(O) = 'Ct^''(0) = as (2.61), therefore (2.72) is 



satisfied. We also require that the regularity condition for and ^'j:'^' at the horizon, 
which yield the ingoing condition for a2'^ and a[^\ 

Fig. 5-10 show Im[(T''^''(aj)] for various wavenumbers g = 2, 5, 10 for each case. 
Let us define the coefficient C{q) as Im[(T^"'^^(a;)] = C{q)/u + 0(1) near uj = 0. 
Then, the coefficient C{q) can be read from the best fitting curves in Fig. 5 - 10 as 
0(2) = 0.41, 0(5) = 0.78, and 0(10) = 1.6 for the case (i), while 0(2) = 0.38, 0(5) = 
0.73, and 0(10) = 1.58 for the case (ii). This implies that the coefficient C{q) 
increases as q increases. 
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Figure 7: (color online) Im[(T(a;)'-"'^^] in the case (i) for q = 10. The best fitting curve is 
cr(i) ~ 0.72 + 1.6/a; 
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Figure 8: (color online) Im[(T(a;)(^''] in the case (ii) for q = 2. The best fitting curve is 
o-(i) ~ 2.7 X 10~^ + 0.38/u> 
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Figure 9: (color online) Im[(7(a;)^^''] in the case (ii) for q = 5. The best fitting curve is 
o-(i) ~ 0.0024 + 0.73/a; 



Finally, we show g-dependence of cr(a;)^^^ in Figs. 11 - 16. a{uj)^'^^ is numerically 



){2) 
F(2) 



obtained by solving Eq. (2.35) under the regularity condition for at the horizon 
and the constant electric field condition. Let us expand the imaginary part of the 
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Figure 10: (color online) Im[(T(a;)*-^^] in the case (ii) for q = 10. The best fitting curve is 
o-(^) ~ 0.49 + 1.58/w. 



Imcr2 




-150 - ■ 



Figure 11: (color online) Im[o"((x')^^^] in the case (i) for q = 2. The best fitting curve is 
C7(2) ~ -1.3 X lO"*^ - 0.010/w 



spatially averaged conductivity cr as 

^ CM _ c»°' + .c<-)w + .^c(^>fa) + ... ^ 

near the origin, w = 0. Then, by definition, we immediately obtain C^^\q) = 0. 
Therefore, the coefficient C{q) is given by 

C(g) = C(°) + e2c(2)(g). (2.74) 

The coefficient C^^\q) can be read from the best fitting curves as in Fig. 11 - 
16 as C(2)(2) = -0.010, C^^\5) = -0.014, and C(2)(io) = -0.020 for the case (i), 
while C(2)(2) = -0.0067, C^^\5) = -0.010, and C(2)(io) = -0.015 for the case 
(ii). Therefore, the lattice effects reduce the coefficient C{q) for any wavenumber q. 
Furthermore, we find that as q increases, the coefficient C^'^\q)^ and therefore C(g), 
decreases. 

By Kramers-Kronig relation, the real part of the conductivity contains a delta 
function peak if the imaginary part of the conductivity contains a pole. All these 
results suggest that, the magnitudes of the zero frequency delta function peak de- 
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Figure 12: (color online) Im[(T((x>)^^)] in the case (i) for g = 5. The best fitting curve is 
£7(2) ~ -2.5 X 10^6 - 0.014/w 
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Figure 13: (color online) Im[(T(cj)'-^^] in the case (i) for q = 10. The best fitting curve is 
cr(2) ~ _3.9 X 10-6 - OmO/u 
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Figure 14: (color online) Im[(T(a;)(^''] in the case (ii) for q = 2. The best fitting curve is 
o-(2) ^ _5,6 X 10"^ - 0.0067/a; 

crease by the lattice effects. This implies that in the holographic superconductor, 
the "superfluid component" of the conductivity decreases by the lattice effects. 

3. Conclusion and Discussion 

We studied the lattice effects on the toy model of holographic superconductor (su- 
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Figure 15: (color online) Im[o"(a;)^^''] in the case (ii) for q = 5. The best fitting curve is 
£7(2) ~ 3.5 X 10"^ - O.OlO/w 
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Figure 16: (color online) Im[(T(u;)''^''] in the case (ii) for q = 10. The best fitting curve is 
c7(2) ~ 0.00015 - 0.015/w. 



perfiuidity) , massive U{1) gauge boson model. Especially we studied how the zero 
frequency delta function peak on the real part of the conductivity is influenced by 
the lattice effects. Our analysis suggests that even though its weight reduces, the 
delta- function peak still remains even after the lattice effects are taken into account. 
This implies that the superfiuid component remains with the lattice. We have seen 
also that, as the wavenumber of the lattice increases, the weight of the delta function 
peak decreases. 

However in order to get conclusive results, clearly we need to study things in 
more great detail. In our toy holographic superconductivity (superfluidity) model, 
we have neglected two important ingredients, the dynamics of the charged scalar 
field and also the gravity. For the charged scalar field, instead of treating it as a 
dynamical field, we have given its VEV by hand as an input. This nonzero VEV 
corresponds to the U{1) symmetry breaking and yields the mass term for the gauge 
boson. Even though we vary this VEV and its radial profile through the several 
parameters of our system, we have seen that the zero frequency delta function peak 
remains. Therefore, we expect that the results will not be modified much even after 
we have taken into account the charged scalar dynamics. However of course, it is 
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better to confirm this point in more explicitly by taking into account the dynamics 
of the charged scalar field [18]. 

We have also neglected the effects of the gravity. There are two important effects 
associated with the gravity dynamics; The first one is the back reaction of the lattice 
effects to the geometry, since we have used the background geometry which does 
not possess the lattice effects. This correction can be calculated perturbatively, as 
is studied, for example, in [15]. This induces the perturbative corrections to the 
background geometry, and how the perturbative correction appears on the geometry 
depends on how we introduce it. Suppose the lattice effects are 0(e), then, depending 
on the lattice effects which appear to the energy-momentum tensor at either 0(e) 
or O(e^), the order of back reaction to the geometry is different. If we introduce the 
perturbative lattice effects on the chemical potential as [15] to the background with 
nonzero chemical potential background, then the back reaction of the lattice effects 
appear as 0{Gm^), then it cannot be neglected unless we take the probe limit. On 
the other hand, if we introduce the perturbative lattice effects by introducing neutral 
scalar field as [16], then its back reaction appears as 0{Gn^^) therefore it can be 
neglected at the leading order in e without taking the probe limit. 

There is also gravity effects on the conductivity calculations. As we quote in the 
introduction, the delta function peak does not appear in the normal phase, where 
U{1) symmetry is preserved. This can be seen for example, by the fact that without 
gravity effect, on the normal phase, the equations of motion for the gauge boson 
admits only trivial constant solution^. Once we take into account the gravity, we 
can see the delta function peak, originated from the translational invariance of the 
system. In [16], by taking into account the gravity effect, it is shown that the zero 
frequency delta function peak becomes flatten once we take into account the lattice 
effects. It is interesting to see how the results in this paper are influenced if we take 
into account the gravity effect. 

Finally even if we take into account all of above effects, it is not clear if the 
perturbative analysis, as we have done in this paper, is enough to give us conclusive 
results. It is possible that there are non-perturbative corrections to the conductivity 
by the lattice effects, which significantly influence the delta function peak. In such 
cases, we may have to rely fully on the numerical analysis. We left these open 
questions for the future projects. 
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Figure 17: (color online) Re[(T(tj)(^)] in the case of di = 5, d2 = 0, q = 2. 
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Figure 18: (color online) Re[(T(Li;)(^)] in the case of di = 5, d2 = 0, q = 5. 
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A. A real part of the conductivity 

In this appendix, we summarize the numerical data for the real part of the conduc- 
tivity. In Figs. 2 and 4, we show the Re[(T(ci;)^°)] in the case of di = 5, d2 = and 
di = 5, d2 = 2, respectively. It is clear that there is an energy gap approximately in 
the region < w < 3 for both cases. 

In Figs. 17 - 22, Re[(T(a;)*^*^] {i = 1, 2) are plotted in the case of di = 5, d2 = and 
di = 5, d2 = 2, respectively for each q = 2,5. In any case, the real part approaches 
a constant in the limit a; — )■ 0. 

Note that even though Re[cr(a;)*^*-'] < (z = 1, 2), since we have small but nonzero 
Re[(T(w)(°)] > 0, as long as we consider the perturbative analysis, the real part of the 
conductivities are always positive. 
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Figure 19: (color online) Re[(T(a;)(^)] in the case of di = 5, d2 = 2, q = 2. 
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Figure 20: (color online) Re[(T(a;)(^)] in the case of di = 5, d2 = 2, q = 5. 
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Figure 21: (color online) Re[(T(Li;)(2)] in the case of di = 5, d2 = 0, q = 2. 
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